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A quasi-can ordeal sequence of finite dimension al quantizations has been 
found which has canonical quantization as its hxnit. In order to demonstrate 
its practical utility and its numerical convergence, this formalism is applied 
to the eigenvalue and "eigenfunction” problem of several harmonic and an- 
harmonic oscillators. 
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I. Introduction 

Harmonic and an harmonic oscillators have long been used to illustrate 
new approximation techniques. Here they are used to demonstrate the appli- 
cation of an approximation procedure based upon the approach of a sequence 
of discrete non-canonical quantizations to the standard canonical qu ant iz ac- 
tion limit . 

Consider a discrete one dimensional space in which the coordinates* of 
allowed positions are interger multiples of a fundamental scale parameter a 
having dimensions of length. Although it is well known that the canonical 
commutation relation 

[Q,p] = a/ (i) 

does not admit finite dimensional matrix representations/ one can ask whether 
a limit procedure exists such that sequences of matrices {Qn,Pn} satisfy in 
some way 

Bn. [<?*,P„] = A/ (2) 

in t he w eak sense. We have found the answer to be in the affirmative 3 and 
therefore briefly sketch the theoretical analysis and apply the formalism to the 
numerical eigenvalue and eigenstate problems of harmonic and an harmonic 
oscillators. 

II. Brief Analysis 

The quantum mechanical scalar product of two wave functions in the 
Schroedinger representation can be written in the form 

(¥,$)= f [ V(q)6(q- q f )$(q')dqdq f , (3) 

a form equivalent to the traditional one. With respect to this scalar product 
the Schroedinger coordinate and momentum representatives have the form 

<?(?> «0 = - q')> (4) 

-P(?> ?') = -ih^6{q ~ ( 5 ) 
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On the other hand in a finite discrete space representation of the canonical 
coordinates and momenta one defines the Schroedinger representative of the 
coordinate operator by the diagonal matrix 


/ (JV - 1)/2 0 . 


Qn = a 


1/2 0 

0 . - 1/2 


0 \ 


> 


( 6 ) 


V o 


0 —(N — l)/2 / 


or equivalently in component notation by 


Qn» ~ ra ^ r ,- 


( 7 ) 


Here a defines the distance between neighboring points in space and N is a 
power of 2. 

The components of the momentum operator are taken from the discrete 
Fourier transform of the diagonal part of Qn, 



E 






00 


Thus Pn is a Toeplitz matrix. The expressions (4) and (5) follow from a 
proper definition of the state space scalar product and limiting procedure. 
The state space scalar product is defined by 


= i X) * k * k ( 9 ) 

a k=-{N-l)/2 


where the bar denotes complex conjugation. With respect to this scalar 
product it is convenient for the sake of the weak limiting procedure to write 
the rth eigenvector of Qn as either 

( 10 ) 

r r 
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or equivenently expressed In terms of the finite discrete Fourier transform as 


1 


XN-\)f2 


t - jf £ 

iT *=— <AT— 1)/2 

Using the latter one can verify that 




t i (W-X)/2 

r ' a r, tfa kss _^_ 1)f2 
One can now define the refinement or weak limit as that in which 


'BrM. 


( 11 ) 


(12) 


JVa — *■ oo, a — ► 0, ra — *• g, sa — g 7 , (13) 

so that the matrix product summations carry over into integrals (in the same 
way that a Fourier series can be carried over into a Fourier integral). Note 
that this approach differs from those of others in which the domain of the 
right hand side of (2) is restricted to a subspace of the Hilbert space. 3 

For a small enough and N, r, s large enough one can get as close as desired 
to any real valued q or q\ In this case using the right hand side of (12) one 
can see that in the refinement or weak limit 

(?>?) — %-? 0- (I 4 ) 

Similarly one can see that (2)and (7) and (8) have the weak limits 

(15) 

(?,(??) - (16) 

(?, p ?> ~ ih T q s{q - ( 17 > 

Note that with all of these weak limits the factor - associated with the 

A 

implied metric tensor in the scalar product (9) is essential. 

The existence of these limi ts is considered sufficient justification for inves- 
tigating the practical utility of this finite dimensional quasi-canonical quan- 
tization. Hence we investigate the eigenvalue problem of several of the osdl- 
lators described by the Hamiltonian 

pa „ Q* o* 

E=t— + FQ + + A (18) 

2m 2 4 


I 

I 

I 
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III. Oscillator Eigenstates and Eigenvalues 


Given (8) and (7) the Hamiltonian (18) can be written in the form 
m 1 ( (2 irhk) 2 . „ . Tnur 1 f x3 , mX , x4 \ \ 0 L (r — 5 ) 

= * + + tW ) “P r 2 **— 

(19) 

It now remains to make proper choices of N and a and to carry out the 
numerical calculations in a manner compatible with a refinement process. 

Define exponents n and / such that 


N = 2 n , 0 = 002 -°, L = ao2 n - a t n-ci> 0, (20) 

where a$ is a length scale appropriate to the problem and where L gives the 
physical size of our one dimensional space. Clearly n sets the dimensionality 
of the matrix and a the refinement. The size of the space is Z, while the rth 
eigenvector 9 corresponds to the physical coordinate roo2 n-a . 


A. The Harmonic Oscillator 


For the case of the harmonic oscillator {F — X = 0) it is convenient to 
choose 

— (-M ■ <»> 


With the choices (19), (20), and (21) the discretized dimensionless form for 
the Hamiltonian is 


H_ 

huf 


£ 

*=— (//— 1)/2 


( (2 irk) 


+ 


2 2 a- 


1 ) eIP { _ 2’ r *(“jv")} M 


( 22 ) 


One has different approximations for different choices of n and a . Fig- 
ure I* illustrates the first 4 normalized ”wave n eigenvectors for the discrete 
Harmonic Oscillator ( F = A= 0) with n = 8, ot = 1. Amazingly, for the first 
two eigenvectors of ^ the absolute value of the error between a component 
of the eigenvector and the exact corresponding eigenfunction solution to the 
Schroedinger equation is less that 10” 3 for the ground state and 10“ 2 for the 
first excited state at any of the allowed positions in the discrete space. 
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PROBABILITY AMPIZTOTE 



FIGURE I. HARMONIC OSC WAVEVECTORS (8,0) 
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Table I illustrates the numerical convergence of the 1st four harmonic 
oscillator eigenvalues as a function of matrix dimensionality N and a. 



N = 2,a = 0 

Table I 
N- 4,a = 0 

H.O. 

N = 8, a = 0 

N= 16, a = 

Eo 

1.3587 

.5410 

.50018 

.5000001 

E i 

1.3587 

1.2186 

1.4961 

1.499984 

Ei 

— 

3.1156 

2.5241 

2.500205 

e 3 

— 

3.7933 

3.3582 

3.497943 


B. The Asymmetric Oscillator 

For the case of the asymmetric oscillator we take the choice 

1 


ao 


-(£)■- 


SF F 

A — 


ma o 


ma$' 


(23) 


The components of the dimensionless Hamiltonian now have the form 

E 


TTf ( / 

j£= £ I 

* 30 k=-(N-lV2 \ 


1 h 2 r — 8 • 2™ r 2 H 

23a-3n+l 1 1 4 


— ^ exp | — 2ick— 


zA 

N 


Figure II illustrates the asymmetric potential while Figure III illustrates 
the solutions to the eigenvector problem for this asymmetric oscillator. Of 
interest is the "trapping” in the virtual potential well occuring with the 4th 
energy level, an effect not easily accounted for with other approximation 
techniques. 

Table II illustrates the numerical convergence of the first four eigenvalues 
using several choices of N and ot. 


Table n ASYM OSC. 



AT = 2, a = 0 AT = 4,a = 0 

N = 8, a = 0 

N = 16, a = 1 

N = 64, a = 2 

Eo 

-.7507 -9.3999 

-19.1162 

-19.6787 

-19.6803 

Ex 

-1.2493 -3.2519 

-16.5510 

-15.6227 

-15.6025 

Ei 

-0.6331 

-9.1040 

-11.7496 

-11.7147 

Eo 

- 2.0159 

-8.9120 

-8.6628 

-8.65986 
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q/ a 


— v(q) 

FIGURE II. ASYMMETRIC OSC POTENTIAL 
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PR OB. AXCPLTTtfDE 



_0 -7 -0 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 

q/a 

— — E0=- 19.11 — El =—18.48 E3 =-9.104 _a ~ E4=-8.912 


FIGURE III. ASYM OSC WAVEVECTORS (8,0) 
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C. The Quartic Oscillator 

In the case of the quartic oscillator (F=w = 0) the choice 


_ ( 2irH V 

[mAW/ 


suffices. In this case the components of the dimensionless Hamiltonian have 
the form 

§ = „X + exp {- 2t ^} 

The quartic oscillator is included here because it has been a frequent subject 
of study and its eigenvalues are numerically established. Table III Illustrates 
the numerical convergence for the firet four eigenvalues and gives the "exact” 
eigenvalues as numerically established by others 466 using distinct techniques. 


Table III QUARTIC OSC 



N = 2,ot = 0 N = 4,a = 0 

N = 8, a = 0 

N = 16, a = 1 

"exac£" 

; 


1.2493 

.4498 

.4217 

.420805 

.420804 

! 

Ex 

1.2493 

1.2760 

1.5587 

1.50790 

1.50790 


Ei 

3.0895 

2.9422 

2.95880 

2.95880 

i 

Et 

3.9157 

3.8709 

4.62128 

4.62122 

* 


IV. Summary. 

The purpose here has been to demonstrate the numerical application 
of the quasi-canonical quantization procedure to 1 dimensional oscillatory 
systems. With this application it is clear that numerical convergence occurs 
for each of the potentials. 

The calculation procedure is computationally straightforward since the 
momentum matrix (8) can be conveniently calculated using well established 
signal processing techniques. 7 In particular, fast Fourier transform (FFT) 8 
techniques allow for rapid determination of "state” vectors and quantum 
numbers. 
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The approximation procedure utilized here is based not upon the trun- 
cation of the wave functions but rather on how well the canonical commutar 
tion relation is approximated. In contrast to usual perturbation theories the 
number of allowed states in each order is thesame as the number of allowed 
positions. The convergence of this discrete quantization to a canonical one 
and detailed derivations of results quoted in the Section II will be discussed 
elsewhere. However it is useful to note two important facts about the pro- 
cedure. First, the use of the discrete fourier transform (and numerical use 
of FFT’s 8 ) does not implies the periodicity of space, rather the quantization 
is carried out over a finite region of space. Secondly, the matrix Q // as a 
N=2 n dimensional matrix insures not only the faster speed of the FFT’s uti- 
lized but also that on the discrete scale that Qjy is an invertible matrix for 
which the eigenvalue zero exists only as a refinement or weak limit. Thus 
one should not expect numerical problems in dealing with coulomb like po- 
tentials. 

Finally, with the refinement limit we note that one can get as close as 
one wishes to canonical quantization with discrete space, suggesting that 
quantum theory cannot readily distinguish between discrete and continuous 
space time. In addition, please note that the combination (9) and (11) fd- 
lows a sequence of matrices to serve as the definition of required generalized 
functions, as in (14,15,16,17), somewhat in analogy to the good functions of 

Lighthill 9 . 10 
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